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Chapter 1 



1 . The vectors x + y + z and -x - y + z are in the directions of two body diagonals of a 
cube. If 0 is the angle between them, their scalar product gives cos 0 = -1/3, whence 
0 = cos '1/3 = 90° + 19° 28' = 109° 28'. 

2. The plane (100) is normal to the x axis. It intercepts the a' axis at 2a' and the c' axis 
at 2c' ; therefore the indices referred to the primitive axes are (101). Similarly, the plane 
(001) will have indices (Oil) when referred to primitive axes. 

3. The central dot of the four is at distance 




cos 60° a 

a = a ctn 60 = — 1 = 

cos 30° V3 



from each of the other three dots, as projected onto the basal plane. If 
the (unprojected) dots are at the center of spheres in contact, then 



a 2 = 



V V3y 



+ 



UJ 



or 




C -f-= 1.633. 
a v 3 
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Chapter 2 



1. The crystal plane with Miller indices hkf is a plane defined by the points ai/h, a 2 /k, and a 3 II . (a) 
Two vectors that lie in the plane may be taken as a,/h - a 2 /k and a, / h — a , / i . But each of these vectors 
gives zero as its scalar product with G = ha, + ka, + fa 3 , so that G must be perpendicular to the plane 
hkf . (b) If n is the unit normal to the plane, the interplanar spacing is n • a,/h . But n = G / I G I , 
whence d(hk£) = G a, / hIGI = 2rc / I Gl . (c) For a simple cubic lattice G = (2nl a)(hx + ky + f, z) , 
whence 

1 _ G 2 _ h 2 +k 2 + 1 2 
a 

— >/3a —a 0 
2 2 

2. (a) Cell volume a, -a 2 xa 3 = -^-v/3a ^-a 0 

0 0 c 



= — y/3 a 2 c. 

2 

x y z 

, u , K O a 2 xa 3 4n 1 ft 1 
(b) b, = 2n — = _ v3a —a 0 

la,-a 2 xa 3 l V 3a"c 2 2 

0 0 c 

2n 1 ^ ^ 

= — (— ^=x + y), and similarly for b 2 , b 3 . 
a V3 




(c) Six vectors in the reciprocal lattice are shown as solid lines. The broken 
lines are the perpendicular bisectors at the midpoints. The inscribed hexagon 
forms the first Brillouin Zone. 



3. By definition of the primitive reciprocal lattice vectors 



Vbz _ {2nf (a 2 x, Ha i xa,)x(a,xa,) _ ^ , , . x , , 

I (a, -a 2 xa 3 )' I 

= (27t) 3 /V c . 



For the vector identity, see G. A. Korn and T. M. Korn, Mathematical handbook for scientists and 
engineers, McGraw-Hill, 1961, p. 147. 

4. (a) This follows by forming 
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l-exp[-iM(a- Ak)] l-exp[iM(a • Ak)] 
1 - exp[-i(a • Ak)] l-exp[i(a • Ak)] 

1 - cos M(a • Ak) _ sin 2 \ M(a • Ak) 
l-cos(a-Ak) sin 2 |(a-Ak) 



(b) The first zero in sin — Ms occurs for e = 2jt/M. That this is the correct consideration follows from 
2 



sin M(7th + — s) = sin 7tMh cos — Ms + cos 7tMh sin — Ms. 

2 ' 2 ^7— 2 



zero, 
as Mh is 
an integer 



5. S (VjVjV,) = f E e 



-27ti(x j v 1 +y j v 2 +z j v 3 ) 



Referred to an fee lattice, the basis of diamond is 000; . Thus in the product 

4 4 4 



S( v iV 2 v 3 ) = S(fcc lattice) x S (basis) , 



we take the lattice structure factor from (48), and for the basis 



s 



(basis) = 1 + e 



. 1 
1 — 

2 



7t( Vi + V 2 + v 3 )- 



Now S(fcc) = 0 only if all indices are even or all indices are odd. If all indices are even the structure factor 
of the basis vanishes unless Vi + v 2 + v 3 = 4n, where n is an integer. For example, for the reflection (222) 
we have S(basis) = 1 + e _l311 = 0, and this reflection is forbidden. 

/• CO . ^ . 

6. f G = 47rr“ (7ra 0 ~ Gr) sin Gr exp (-2r/a 0 ) dr 

= (4/ G 3 a 0 3 )J dx x sin x exp (-2x/Ga 0 ) 

= (4/G 3 a 0 3 ) (4/Ga 0 )/(l + r/G 2 a 0 2 ) 2 
16/(4 + G 2 a 0 2 ) 2 . 

The integral is not difficult; it is given as Dwight 860.81. Observe that f = 1 for G = 0 and f oc 1/G 4 for 

Ga 0 » 1. 

7. (a) The basis has one atom A at the origin and one atom B at —a. The single Laue equation 

2 

a • Ak = 2 n x (integer) defines a set of parallel planes in Fourier space. Intersections with a sphere are 
a set of circles, so that the diffracted beams lie on a set of cones, (b) S(n) = fA + fB e _llm . For n odd. S = fA - 
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f B ; for n even, S = f A + f B . (c) If fA = f B the atoms diffract identically, as if the primitive translation vector 
were — a and the diffraction condition (—a • Ak) = 2 71 x (integer). 
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Chapter 3 



1. E = (M 2 /2M) (In/X) 2 = (M 2 /2M) (ti/L ) 2 , with X = 2L. 

2. bcc: U(R) = 2Ns[9.1 14(a/R) u - 12.253(a/R) 6 ]. At equilibrium R 0 6 = 1.488a 6 , and 
U(R 0 ) = 2Ne(- 2.816). 

Ik: U(R) = 2Ns[12.132(a/R) 12 - 14.454(a/R) 6 ]. At equilibrium R 0 6 = 1.679a 6 , and 

U(R 0 ) = 2Ns( — 4.305). Thus the cohesive energy ratio bcc/fcc = 0.956, so that the fee structure is 
more stable than the bcc. 



3. | U | = 8.60 Ns 

= (8.60)(6.02 x 10 23 ) (50 x 10 16 ) = 25.9 x 10 9 erg/mol 
= 2.59 kJ/mol. 



This will be decreased significantly by quantum corrections, so that it is quite reasonable to find the same 
melting points for H 2 and Ne. 

4. We have Na — > Na + + e - 5.14 eV; Na + e — > Na~ + 0.78 eV. The Madelung energy in the NaCl 
structure, with Na + at the Na + sites and Na~ at the Cl sites, is 



10x2 



ae _ (1.75) (4.80 x 10 1U ) 



= 11.0 x 10 12 erg, 



R 



3.66 x 10- 



or 6.89 eV. Here R is taken as the value for metallic Na. The total cohesive energy of a Na 1 Na~ pair in the 
hypothetical crystal is 2.52 eV referred to two separated Na atoms, or 1.26 eV per atom. This is larger than 
the observed cohesive energy 1.13 eV of the metal. We have neglected the repulsive energy of the Na + Na~ 
structure, and this must be significant in reducing the cohesion of the hypothetical crystal. 

5a. 



U(R) = N 



( A aq^ 3 

y R" R^ j 



; a = 2 log 2 = Madelung const. 



In equilibrium 



au 

dR 



= N 



nA aq 



2 3 



+ ■ 



V R o" R„ 2 y 



= 0 ; 

aq 



and 



U(R 0 ) = -^31(1-4). 

R 0 n 
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b. U(R 0 -R,8) = U (R„ ) + ||^-R 0 (R,8) 2 + . . . , 
bearing in mind that in equilibrium (0U/0 R) r = 0. 



f d 2 \ J A 



dR 2 



= N 



^n(n + l)A 2aq 2 ^ 



R„ 



v R 0 



n+2 



R 



= N 



o J 



(n + l)aq 2 2aq 2 



R n 



R 



o J 



For a unit length 2NR 0 = 1, whence 



f d 2 \ 



V^R 2 7 



C = R 0 2 ®A 



R 



2R 



o 



dR 2 



(n-l)q 2 log 2 

R « 2 



6. For KC1, X = 0.34 x 10 8 ergs and p = 0.326 x 10 8 A. For the imagined modification of KC1 with the 
ZnS structure, z = 4 and a = 1.638. Then from Eq. (23) with x = R 0 /p we have 

x 2 e~ x = 8.53 x 1(T 3 . 



By trial and error we find X — 9.2, or R 0 = 3.00 A. The actual KC1 structure has R 0 (exp) = 3.15 A . For 
the imagined structure the cohesive energy is 



u=-' aq " 



R„ 



_P 

R 



or -^=-0.489 



o J 



in units with R 0 in A. For the actual KC1 structure, using the data of Table 7, we calculate — — = —0.495, 

q" 

units as above. This is about 0.1% lower than calculated for the cubic ZnS structure. It is noteworthy that 
the difference is so slight. 

7. The Madelung energy of Ba + CT is -ae 2 /R 0 per ion pair, or -14.61 x 10~ 12 erg = -9.12 eV, as compared 
with -4(9.12) = -36.48 eV for Ba ++ O . To form Ba + and CT from Ba and O requires 5.19 - 1.5 = 3.7 eV; 
to form Ba ++ and O'" requires 5.19 + 9.96 - 1.5 + 9.0 = 22.65 eV. Thus at the specified value of Ro the 
binding of Ba + CT is 5.42 eV and the binding of Ba 1 1 O " is 13.83 eV; the latter is indeed the stable form. 

8. From (37) we have exx = SnX x , because all other stress components are zero. By (51), 

3S, 1 =2/(C m -C, ! )+1/(C„ + C, ! ). 

Thus Y — (C, +C 12 C„ -2C, ! ! )/(C„ +C,,); 

further, also from (37), eyy = S 2 iX x , 

whence a^e^/e^ = S 21 /S n = -C 12 /(C n +C 12 ). 

9. For a longitudinal phonon with K || [1 1 1], u = v = w. 
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co 2 p = [Q , + 2C 44 + 2(C 12 + C 44 )] K 2 /3 , 
or v = co/K = [(C n + 2C 12 + 4C 44 /3p)] 1/2 

This dispersion relation follows from (57a). 

10. We take u = - w; v = 0. This displacement is _L to the [111] direction. Shear waves are degenerate in 
this direction. Use (57a). 

11. Let e xx =-e yy =X e in (43). Then 



U = XC n (%e 2 +Xe 2 )-XC 12 e 2 

= K[K(C 11 -C 12 )]e 2 



so that 



f 5 2 U^ 


= N 


Ur 2 ; 


R 0 



n(n + 1)A 



R 



n+2 



constant. 



2aq 2 ^ 



(n + l)aq 2 

V Ro 3 



2aq 2 2 

~R7~, 



is the effective shear 



12a. We rewrite the element ay = p - °y(L + p - q) as ay = p - X’ 8 y, where X' = X + p - q, and °y is the 
Kronecker delta function. With X' the matrix is in the “standard” form. The root X' = Rp gives X = (R - l)p 
+ q, and the R - 1 roots X' = 0 give X = q - p. 

b. Set 



u (r, t) = u 0 e i[(KM)(x+y+z) -“ t] ; 

v(r,t) = v 0 e l[ J ; 

w(r,t) = w 0 e l[ ] , 

as the displacements for waves in the [111] direction. On substitution in (57) we obtain the desired 
equation. Then, by (a), one root is 

co 2 p = 2p + q = K 2 (C n + 2C 12 + 4C 44 )/ 3, 

and the other two roots (shear waves) are 

co 2 p = K 2 (C 11 -C 12 +C 44 )/3. 

13. Set u(r,t) = u 0 e l<K - r ~ t) and similarly for v and w. Then (57a) becomes 

co 2 pu 0 =[C 11 K y 2 +C 44 (K y 2 +K z 2 )]u 0 
+ (C 12 +C 44 ) (K x K y v 0 +K x K z w 0 ) 

and similarly for (57b), (57c). The elements of the determinantal equation are 
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M 11 =C 11 K x 2 +C 44 (K y 2 +K z 2 )-«V; 

M 12 = (C 12 + C 44 )K x K y ; 

M 13 =(C 12 +C 44 )K x K z . 

and so on with appropriate permutations of the axes. The sum of the three roots of 0) p is equal to the 
sum of the diagonal elements of the matrix, which is 

(C M +2C 44 )K 2 , where 

K 2 =K x 2 +K y 2 +K z 2 , whence 
Vi 2 +V 2 2 +V 3 2 =(Cn+2C 44 )/p, 

for the sum of the (velocities) 2 of the 3 elastic modes in any direction of K. 

14. The criterion for stability of a cubic crystal is that all the principal minors of the quadratic form be 
positive. The matrix is: 

C„ C 12 C 12 

C12 Ci i C 12 

Ci 2 C i2 C„ 



The principal minors are the minors along the diagonal. The first three minors from the bottom are C44, 
C44 2 , C44 3 ; thus one criterion of stability is C44 > 0. The next minor is 

Cn C 44 3 , or C„ > 0. Next: C 44 3 (C„ 2 - C 12 2 ), whence |C 12 | < C„. Finally, (C„ + 2C 12 ) (C„ - C 12 ) 2 > 0, so 
that Cn + 2C [2 > 0 for stability. 
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Chapter 4 



1 2 

la. The kinetic energy is the sum of the individual kinetic energies each of the form — Mu s . The force 

between atoms s and s+1 is -C(u s - u s+ i); the potential energy associated with the stretching of this bond is 

1 2 

— C(U s — U J+1 ) , and we sum over all bonds to obtain the total potential energy. 



1 2 * 1 . . 22 

b. The time average of — Mu s IS — Mco“ir. In the potential energy we have 



u s+1 = ucos[cot-(s + l)Ka] = u{cos(cot-sKa)-cos Ka 
+ sin (cot - sKa) • sin Ka } . 

Then u s -u s+1 =u {cos(cot-sKa)-(l-cos Ka) 

- sin (cot - sKa) • sin Ka} . 

We square and use the mean values over time: 

< cos 2 > = < sin 2 > = — ; < cos sin > = 0. 



Thus the square of u{ } above is 

^-u 2 [1-2cos Ka + cos 2 Ka + sin 2 Ka] = u 2 (l-cos Ka). 

The potential energy per bond is — Cu 2 (l — COS Ka), and by the dispersion relation to 2 = (2C/M) (1 - 



1 



2 

2 2 



cos Ka) this is equal to — Mco U . Just as for a simple harmonic oscillator, the time average potential 
energy is equal to the time-average kinetic energy. 

2. We expand in a Taylor series 



+ 





+— p 2 a 2 


( a 2 u "j 


y 8XJ 


s 2 P 





u(s + p) = u(s) + pa 

On substitution in the equation of motion (16a) we have 

■ = ( z pVc n ) 



M 52u _rv „2„2^ .5 2 U 



at 2 p>o r p/ a x 2 

which is of the form of the continuum elastic wave equation with 
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v 2 = IVT 1 2 p 2 a 2 C . 



3. From Eq. (20) evaluated at K = jt/a, the zone boundary, we have 

-orMjit = -2Cu ; 

-co 2 M 2 v = -2Cv . 

Thus the two lattices are decoupled from one another; each moves independently. At ro 2 = 2C/M 2 the 
motion is in the lattice described by the displacement v; at co 2 = 2C/Mi the u lattice moves. 



4. 


co 2 


2 


A 2 


sin pk 0 a 

— (1-cos pKa) 






M 


p>0 


pa 


5co 2 


2A 




sin pk 0 


a sin pKa 






2 


GK " 


" M 


p>0 



— (cos (k 0 - K) pa - cos (k 0 + K) pa) 

When K = k 0 . 



3co 2 

8K 



— 2 (1-cos 2k 0 pa) , 

M P>0 



which in general will diverge because 2 1 — > oo. 

P 



5. By analogy with Eq. (18), 



Md 2 u s /dt 2 = C,(v s -u s ) + C 2 (v s _, -u s ); 
Md 2 v s /dt 2 =C,(u s -v s ) + C,(u s+] - v ). whence 
-co 2 Mu = C 1 (v-u) + C 2 (ve~ lKa -u); 

-co 2 Mv = C 1 (u-v) + C 2 (ue lKa -v) , and 



(C, +C 2 )-Mco 2 -(Cj + C 2 e~ iKa ) 
-(Cj +C 2 e iKa ) (C, +C 2 )-Mco 2 






(20C/m;' /i For Ka = 0, co 2 =0and 2(C 1 +C 2 )/M. 

For Ka = n, co 2 = 2Q/M and 2C 2 /M. 




(2C/1*) 



6. (a) The Coulomb force on an ion displaced a 
distance r from the center of a sphere of static or rigid conduction electron sea is - e 2 n(r)/r 2 , where the 
number of electrons within a sphere of radius r is (3/4 jtR 3 ) (4rcr 3 /3). Thus the force is -e 2 r/R 2 , and the 
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• 2 3 1/2 2 3 1/2 

force constant is eVR . The oscillation frequency o>d is (force constant/mass) , or (eVMR ) . (b) For 

sodium M —4x10 “ 3 g and R — 2x10 8 cm; thus C0 D — (5 x 10 10 ) (3 x 10 46 )*/ 2 
— 3 x 10 13 S 1 (c) The maximum phonon wavevector is of the order of 10 s cm -1 . If we suppose that ro 0 is 

associated with this maximum wavevector, the velocity defined by a>o/K max « 3 x 10 5 cm s _1 , generally a 
reasonable order of magnitude. 

7. The result (a) is the force of a dipole e p u p on a dipole e 0 u 0 at a distance pa. Eq. (16a) 
becomesco 2 = (2/ M)[y(l- cos Ka) + Z (-l) p (2e 2 /p 3 a 3 )(l-cos pKa)] . 



At the zone boundary oo 2 = 0 if 



1 + a Z (-l) p [l-(-l) p ]p 3 =0 , 



or if a Z[1 — ( — l) p ]p 3 = 1 . The summation is 2(1 + 3“ 3 + 5~ 3 + ...) = 2.104 and this, by the properties of 
the zeta function, is also 7 C, (3)/4. The sign of the square of the speed of sound in the limit Ka « 1 is 
given by the sign of 1 = 2a Z (— l) p p 3 p“, which is zero when 1 - 2 _1 + 3 _1 - 4 _I + . . . = l/2o. The series 

p>0 

is just that for log 2, whence the root is a - 1/(2 log 2) = 0.7213. 



4-3 




Chapter 5 



1. (a) The dispersion relation is CO = C0 m I sin — Ka! . We solve this for K to obtain 

K = (2/a) sin '(co/co^J, whence dK/dco = (2/a)(co m ~ — C0“) 17 ~ and. from (15), D(co) 
= (2L/7ta)(cO m ~ — CO 2 ) 1/_ . This is singular at ro = ro m . (b) The volume of a sphere of radius K in 
Fourier space is Q = 4nK' / 3 = (4tc / 3)[(co 0 — Co)/ A] 3/2 , and the density of orbitals near (B 0 is 

D(co)= ( L/ 2 TC ) 1 I dQ/dco 1= (L/2 tt) 3 (2tt / A 3/ ~)(co 0 — C0) 1/_ , provided to < ro 0 . It is apparent that 
D(cb) vanishes for to above the minimum coo- 



1 2 3 1 1 

2. The potential energy associated with the dilation is — B(AV/V)”a~ « — k B T . This is — k B T and not 
3 

— k B T , because the other degrees of freedom are to be associated with shear distortions of the lattice cell. 

Thus <(AV) 2 >=1.5xlO“ 47 ;(AV) rms = 4.7x1 (T 24 cm 3 ; and (AV) rms /V = 0.125 . Now 

3Aa/a ~ AV/V , whence (Aa) rms / a = 0.04 . 



3. (a) < R >= (K/2pV)Zco 1 , where from (20) for a Debye spectrum Eco 1 

= 1 dco D(co)co ' = 3Vco d " / 47T 3 V 3 , whence < R~ > = 3Kco D ~ / 87 C pv 3 . (b) In one dimension from 
(15) we have D(co) = L/7CV , whence jdco D(co) CO 1 diverges at the lower limit. The mean square 

strain in one dimension is < (SR/5x)“ >= — EK~U 0 2 = (M/2MNv)EK 

= (M/2MNv) (K d 2 / 2) = Mco d 2 / 4MNv 3 . 

4. (a) The motion is constrained to each layer and is therefore essentially two-dimensional. Consider one 
plane of area A. There is one allowed value of K per area (2 ji/L) 2 in K space, or (L/2 ji) 2 = A/4jt 2 allowed 
values of K per unit area of K space. The total number of modes with wavevector less than K is, with co = 
vK, 

N = ( A/4tT ) (ttK 2 ) = Aco 2 / 4tcv 2 . 

The density of modes of each polarization type is D(ro) = dN/dw = Aco/2tiv 2 . The thermal average phonon 
energy for the two polarization types is, for each layer, 

TT , -f«>D Aco /zco 

U = 2 D(co) n(co,T) /?co dco = 2 v dco, 

Jo Jo 2 ttv 2 exp(hco/T)-l 

j* 

where ca D is defined by N = D(co) dco . In the regime /zoo D » T , we have 

J D 



U = 



2At 3 



2izv 2 fi 2 Jo e x -l 



dx. 



5-1 




Thus the heat capacity C = k B SU/St oc T” . 

(b) If the layers are weakly bound together, the system behaves as a linear structure with each plane as a 
vibrating unit. By induction from the results for 2 and 3 dimensions, we expect C oc T . But this only 
holds at extremely low temperatures such that X « /tCO D ~ /zvN layer / L , where Ni ayer /L is the number of 
layers per unit length. 

5. (a) From the Planck distribution < n > +— = — (e x + 1) /(e* — 1) = — COth (x/2) , where 

X = Mco/k B T . The partition function Z = e x/ ~ X C = C x/_ /(l — e x ) = [2sinh (x/2)] 1 and the 
free energy is F = k B T log Z = k B T log[2 sinh(x/2)]. (b) With ro(A) = o>(0) (1 - yA), the condition 

oF/GA = 0 becomes BA = yZ Hot coth (Mco/2k B T) on direct differentiation. The energy 

< n > Mco is just the term to the right of the summation symbol, so that BA = yU (T) . (c) By definition 
of y, we have 8co / CO = — y5V/V , or d log CO = — 8 d log V . But 0 oc C0 D , whence 
d log 0 = -y d log V . 
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Chapter 6 



K 2 , 2 

1. The energy eigenvalues are 8 k = k . The mean value over the volume of a sphere in k space is 



K 2 /k 2 dk-k 2 3 M 2 , 2 3 

< s > = = k ~ = — s F . 

2m /k"dk 5 2m 5 



The total energy of N electrons is 



U 0 =N--8 F . 



2a. In general p = -5U/0V at constant entropy. At absolute zero all processes are at constant entropy (the 

, TT TT 3 xt 3 M 2 f 3n 2 N V /3 

Third Law), so that p = — dU 0 /dV, where U 0 =— Ns F = — N , whence 

5 5 2m V , 



p = . (b) Bulk modulus 

3 V 

B = -v^.=vf-2Lo + A^L2.To + f2YT» = 12To 

dV 3V ! 3V dvj 3 V 1.3J V 9 V' 

(c) For Li, 



^ = |(4.7 x 10 22 cnT 3 )(4.7 eV)(1.6 xlO 12 erg/eV) 

= 2. 1 x 10 1 1 erg cm 3 = 2. 1 x 10 1 1 dyne cm 2 , 
whence B = 2.3 x 10 11 dyne cm -2 . By experiment (Table 3.3), B = 1.2 x 10 u dyne cm -2 . 

rco 1 

3. The number of electrons is, per unit volume, n = I ds D(s) • ^ , where D(e) is the density 

^ e T 1 



of orbitals. In two dimensions 



m f°°^ 1 

^Jo £ Q te-v)h _|_ i 



-jCp + t log (1 + e^ /T )), 



where the definite integral is evaluated with the help of Dwight [569.1]. 



2xl ° 33 m57 

4a. In the sun there are — — 10 nucleons, and roughly an equal number of electrons. In a 

1.7 xlO' 24 

white dwarf star of volume 
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— (2xl0 9 ) 3 *3xl0 28 cm 3 



the electron concentration is 



10 57 

3 x 10 28 



« 3 x 10 28 



cm 3 . Thus 



£ f 



VL 2 i i 

(37l 2 n) 2 ' 3 * —10 " 7 -10 J) » —10 7 ergs, or * 3.10 4 eV. (b) The value of k F is not 

2m 2 2 



affected by relativity and is « n l/3 , where n is the electron concentration. Thus S F — Mck p — Me 3 Vn. (c) A 

6 183 • 38 

change of radius to 10 km =10 cm makes the volume « 4 x 10 cm and the concentration ~ 3 x 1 0 cm 
3 . Thus S F *10 - 7 (3.10 IO )(10 13 ) * 2.10 4 erg * 10 s eV. (The energy is relativistic.) 



5. The number of moles per cm 3 is 81 x 10 3 /3 = 27 x 10 3 , so that the concentration is 16 x 10 21 atoms cm 
3 . The mass of an atom of He 3 is (3.017) (1.661) x 10“ 24 = 5.01 x 10“ 24 g. Thus 
S F =[(1.1 x 1 0 54 ) /l 0 23 ] [(30)(1 6) x 10 21 ] 2/3 *7 x 10 16 erg,orT F *5K. 

6. Let E, v vary as e“ lwt . Then 



V 



eE/m exE 1 + icox 
-ico + (l/x) m 1 + (cox) 2 



and the electric current density is 



j 



= n(-e)v = 



ne 2 x 

m 



1 + icox 
1 + (cox) 2 



7. (a) From the drift velocity equation 

icov x = (e/ m)E x + co c v y ; icov y = (e/ m)E y - co c v x . 



We solve for v x , v y to find 

(co c 2 -co 2 )v x =ico(e/m)E x +co c (e/m)E y ; 

(co c 2 -co 2 )v y =ico(e/m)E y +co c (e/m)E x . 

We neglect the terms in ro c 2 . Because j = n(-e)v = ctE, the components of a come out directly, (b) From the 
electromagnetic wave equation 

c 2 V 2 E = sd 2 E/5t 2 , 

we have, for solutions of the form e l(kz " GA \ the determinantal equation 

2 2i 2 2 

s xx co-ck s xy co = 

s yx co 2 e^-cV 
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Here S xx = 8 yy =1 — CO p ~/co“ and 8 xy = — S yx = i (D c (Q p 2 / CO 3 . The determinantal equation gives the 
dispersion relation. 

8. The energy of interaction with the ion is 

e£" (p/r)47t:r 2 dr = -3e 2 /2r 0 , 

where the electron charge density is -e(3/47tr 0 3 ). (b) The electron self-energy is 

p 2 | o ° dr(4Tcr 3 /3)(47tr 2 )r _I =3e 2 /5r 0 . 

The average Fermi energy per electron is 3ep/5, from Problem 6.1; because N/V = 3 / 4 711',/ , the average 
is 3 (9tc/ 4) K 2 /l 0 mr (| ~ . The sum of the Coulomb and kinetic contributions is 

TT 1.80 2.21 

U = + — 

r s r s 

which is a minimum at 

= or r s =4.42/1.80 = 2.45. 
r r 

s s 

The binding energy at this value of r s is less than 1 Ry; therefore separated H atoms are more stable. 

9. From the magnetoconductivity matrix we have 

c CO c T 

Jy =a yx E x=— T^oE x • 

1 + (co c t) 

For g) c t » 1. we havea yx = a 0 /co c T = (ne 2 T/m )(mc/cBx) = neB/c . 

10. For a monatomic metal sheet one atom in thickness, n ~ 1/d 3 , so that 

R sq « mv F /nd 2 e 2 « mv F d/e 2 . 

If the election wavelength is d, then mv F d ~ H by the de Broglie relation and 

R sq * M/e 2 = 137/c 

in Gaussian units. Now 
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R sq (ohms) = 10 9 c R sq (gaussian) 
« (30) (137) ohms 
*4.1kQ . 
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Chapter 7 



la. The wavevector at the corner is longer than the wavevector at the midpoint of a side by the factor V2. 
As e oc k 2 for a free electron, the energy is higher by (V 2) 2 = 2. b. In three dimensions the energy at a 
corner is higher by (V 3) 2 than at the midpoint of a face. c. Unless the band gap at the midpoint of a face is 
larger than the kinetic energy difference between this point and a corner, the electrons will spill over into 
the second zone in preference to filling up the corner states in the first zone. Divalent elements under these 
conditions will be metals and not insulators. 



2. 8 = M k /2m, where the free electron wavevector k may be written as the sum of a vector K in the 
reduced zone and of a reciprocal lattice vector G. We are interested in K along the [111] direction: from 

Chap. 2, K = (2 n / a) (1,1,1) u, with 0 < U < — , will lie in the reduced zone. 

The G's of the reciprocal lattice are given by 

G = (27t/a)[(h-k + ^)x + (h + k-£)y+ (-h + k + £)z], where 




h,k. 



are 



any integers. Then 8 = 

(27t/a)- [(u + h- k + f)"+(u + h + k-£)“+(u-h + k + (')“].We now 
have to consider all combinations of indices h,k ,£ for which the term in brackets is smaller than 
6[3(l/2) 2 ] or 9/2. These indices are (000); (l 1 l); (lOO), (OK)), and (00l); (100), (010), and 
(001); (111); (HO), (TOT), and (Oil); (110), (101), and (011). 



3. (a) At k = 0 the determinantal equation is (P/Ka) sin Ka + cos Ka = 1. In the limit of small positive P this 
equation will have a solution only when Ka <§C 1. Expand the sine and cosine to obtain in lowest order 







The 



energy 



is 



8 = 



M K / 2m — 1) P / ma . (b) At k = n/a the determinantal equation is (P/Ka) sin Ka + cos Ka = -1. In 
the same limit this equation has solutions Ka = n + 5, where § <gC 1 . We expand to obtain 

(P/7t:)(— 8) + j^— 1 + — 5" = — 1, which has the solution § = 0 and 5 = 2 P/ji. The energy gap is 
E g = (K 2 /2ma 2 )(2tc 5) - (H 2 /2ma 2 ) (4P) . 



4. (a) There are two atoms in the basis, and we label them a and b. Then the crystal potential may be 

f 1 1 1 ' 

x H — a, y + — a, z + — a 
4 4 4 . 



written as U = Uj + U, = U, (r) + U, 
components U G = U 1G + U-, G = U 1G 



and the Fourier transform has 



^ i(G x +G y +G z )/ a ^ 

1 + e y u 



V 



. If G = 2 Ax, then the exponential is 



J 



i— Aa 

e 2 = e 111 = — 1, and U G= , A =0, so that this Fourier component vanishes. Note that the quantity in 



parentheses above is just the structure factor of the basis, (b) This follows directly from (44) with U set 
equal to zero. In a higher order approximation we would go back to Eq. (31) where any non-vanishing 
U G enters. 



7-1 




K 2 

5. Let k = K + i H', A,,, = 

11 2m 



(1 ) 

-G 

U J 



±iG H-H- 



The secular equation (46) is now 



X,-e U 
U X_!-e 



= 0, 



M 2 

and for H < < G we have, with <7 = 8 

2m 



(1 ) 
-G 

U J 



A 



a + iG// • 



h 



2 A 



V 



a’ - 



2m 



C 



2m 



a-iG H 



Vi 



2 \ 



V 



2m 



= U 2 - 



G H 



= U, 



2m 



H 2 = 



Uj--a 

H 2 ^ 2 ' 



2m 



G 



6. U(x,y) = - U[e' <2n/a) (x+y) + other sign combinations of + x ± y]. The potential energy contains the four 
reciprocal lattice vectors (2 jt/a) (+1; +1). At the zone corner the wave function e 1(Il/a) (x+y) is mixed with e _1 

(7t/a) (x+y) rp^ e centra j e q Ua ti ons are 



(X-s)C 



n 

a 



n 

a 



UC 



n _ 
a 



n 

a 



= 0 ; 



(X-s)C 



71 


71 


-UC 


r k 


71 1 






a 


a _ 




La 


aj 



where X = 2^M“/2m^7t/a) . The gap is 2U. 
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Chapter 8 



IT) ^ ] 

la. E d =13.60 eVx x — -6.3x10^ eV 

m s’ 

b. r = a„ x 8x — - 6xl0~ 6 cm 

m* 

c. Overlap will be significant at a concentration 



N = — !— « 10 15 atoms cm 

4n r 3 

3 1 



2a. From Eq. (53), n=(n 0 N d ) 1/ "e Et|/ k ® T , in an approximation not too good for the present example. 



n 0 =2 



( m*k B T i 
27th 2 



3/2 



4xl0 13 cm 3 ; 



—^— = 1.45 ; e~ L45 = 0.23 . 

2k B T 

n = 0.46xl0 13 electrons cm 3 . 

b. R H = — — =-1.3xl0~ 14 CGS units 
nec 

3. The electron contribution to the transverse current is 



j y (e) -nep e 



[i e B 



E. +E, 



for the holes j y (h) = nep h 



' -H„B 

v C 



E. +E. 



A 

J 



Here we have used 



C0„,T = 



fk B 



for electrons; 



Vh = 



pt h B 



for holes. 



c c 

The total transverse (y-direction) current is 

0 = (nep e 2 - pe|T h 2 )(B/c)E x + (nep e + pep h )E y , (*) 

and to the same order the total current in the x-direction is 

j x =(peii h +nep e )E x . 



Because (*) gives 




E y =E x B 



PtV-ntV 1 
PPh+ n Pe c 



we have for the Hall constant 

R _ E y _ 1 PPh 2 ~ n P e 2 
j x B ec (pp h + np e ) 2 ' 

4. The velocity components are V x = Mk x / m t ; V = Hk y /m t ; V z = Mk z /m, ; . The equation of motion 
in k space is M dk/dt = — (e/c) V x B . Let B lie parallel to the k x axis; then 
dk x / dt = 0; dk y / dt = -co / ,k z ; co ( , = eB/m^c; dk z / dt = co t k y ; co t = eB/m t c . We differentiate 
with respect to time to obtain d”k y / dt" = —0), dk z / dt ; on substitution for dk z /dt we have 
d"k y / dt" + (0,0), k y = 0, the equation of motion of a simple harmonic oscillator of natural frequency 

co 0 = (co^cOj) 172 =eB/(m / ,m t ) 1/2 c. 

5. Define Q e = eBx e / m e c; Q h = eBx h / m h c . In the strong field limit Q » 1 the 
magnetoconductivity tensor (6.64) reduces to 



fQe 2 


-q; 1 


o' 


2 


r Qh 2 


Qh 1 


o 2 


q; 1 


q; 2 


0 


+ Pe ^h 


-Qh 1 


Qh 2 


0 


0 

V 


0 


1 


n 


l 0 


0 


1 



We can write nec Q c /B for ne”X e / m e and pec Qh/B for pe X h / m h . The strong field limit for o y 
follows directly. The Hall field is obtained when we set 



j =0 = — (n-p)E x + — + ^~ E 

H l Qe Qh J y 



The current density in the x direction is 



ec n , P c r he 

Jx = — — + ~ E -(n-p)E ; 

B l Q e Q h J P y 



using the Hall field for the standard geometry, we have 



ec n p] (n-p)“ 

Jx=— + — +7 — 

H (Q, Q,J + 

V Qe Qh 
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Chapter 9 

1 . 





b. 



N = 2x F 0 
(27T/k)" 

n = N/L 2 = kp / 27T 
k F = V27tn 

n = — xl0 16 els/cm 2 
8 

k F =0.89xl0 8 cm 1 




banal 




Z# b oH 



1 

3a. In the hep structure there is one atom whose z coordinate is 0 and one at — C . The structure factor of 

2 

„ 2n „ 

the basis for Cj c = Z is 

C 



S G (basis) = 1 + e~ i7t =1-1 = 0, 

Sr c 

so that by the same argument as in Problem 9.4 the corresponding component UG C of the crystal potential 
is zero. 

b. But for U 2G the structure factor is 



S 2Gc (basis) = 1 + e~ i2n = 2. 

c. The two valence electrons can just fill the first BZ. All we need is an adequate energy gap at the zone 
boundary and for simple hex. there is no reason against a gap. 

d. In hep there will be no gap (at least in lowest order) on the top and bottom faces of the BZ, by the 
argument of part a. 
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b. The electron moves in a direction normal to the Fermi surface — more or less in a straight line if the 
Fermi surface is close to planar in the region of interest. The magnetic field puts a wiggle on the motion, 
but the field does not make the electron move in a helix, contrary to the behavior of a free electron. 

6a. 




«»0 



€*-Uo 



Region I: 



M 2 



2m dx 



2 -u 0 



y/ = ey/ 



M 2 k 2 

y/ = A cos kx ; s = U 0 (*) 

2m 



Region II: 



Ir 

2m 



d 2 

dx 2 



y/ = &y/ 



y/ = B e qx ; s = 



mV 

2m 



1 dy/ 

Boundary condition continuous. 

y/ dx 



k tan (ka / 2) = q , 



(*) 






with k and q related to e as above. 

b. The lazy way here is to show that the e’s in the equations marked (*) above are equal when k and q are 
connected by (**), with s = -0.45 as read off Fig. 20. This is indeed so. 

. , 1 . 27te , „ 

7a. A ( — ) = , where S = Jtk F , with k F = 0.75 x 10 cm from Table 6.1, for potassium. Thus 

H McS 



A(— ) - 0.55 x 10 8 G 1 . 

H 137 k F “ e 
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b. 



® c R = v f ; 



_ v F mc _ ftk F c 
eB eB 
- 0.5xl(T 3 cm 



ttR 2 = 0.7x10 6 cm 2 . 



8. Write (17) as H = H () + H l , where H } = (ft / m) k • p . Then (18) is an eigenfunction of H () with 
the eigenvalue 8 n (0) + ft" k" / 2m . In this representation the diagonal matrix element of H { is equal to 
(ft / m) J dV U 0 (r)k-pU 0 (r). In a cubic crystal U () (r) will be even or odd with respect to the 
inversion operation r -> — r , but p is an odd operator. It follows that the diagonal matrix element 

vanishes, and there is no first-order correction to the energy. The function U k (r) to first order in H ] is 



U k (r) = u 0 (r) 



, < jO I H l I nO > 

j e n (0)- ej (0) ’ 



and the energy to second order is 



, , I<n0lk-pl j0>l 2 

s n (k) = 8 n (0) + (ftk) 2 / 2m + (ft / m) 2 Z' — 

j s n (0)-Sj(0) 



The effective mass ratio is the coefficient of ft“ k" / 2m 



or 



m 

m 




m 



knOlpl jO >l 2 

SnW-S^O) ' 



9a. 

JdV w*(r-r n ) w (r-r m ) 

= r'n e lk ' r " e" kr ” f dV y/* (r)^ k (r) 

k k' J k’ 

= N _l I e ik ' (r ”“ rm) 

k 

where the summation is zero unless n = m, when it is equal to N. 
b. w(x-X n ) = N~ l/2 U 0 (x)2:e lk<x x ” \ The summation is 

k 

equal to 
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£ ^i27ip(x-x n )/Na 

P 



1 



-N/2 



i2jip(x-x n )/Na dp 

c 



el ,(x-x n )/a _ e -i,(x-x n )/a _ ^ [^(x-xj/a] 

i27r(x-x n )/Na 7t(x-x n )/Na 



whence 



w(x-x n ) = N 1/2 u 0 (x) 



sin [7i(x-x n )/a] 

rc(X-X n )/a 



10a. j y = Go (Q 1 E x + sE y ) = 0 in the Hall geometry, whence E y = - E x /sQ. 
b. We have j x = ct 0 (Q -2 E x - Q" 1 E y ), and with our result for E y it follows that 

jx=<7 0 (Q~ 2 + s-'Q- 2 )E x , 



whence p = E x /j x = (Q 2 /a 0 ) 



s 

s+r 
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Chapter 10 



d/B 1 

la. — - = — B; this is the London equation. The proposed solution is seen directly to 
dx“ X 



satisfy this and to satisfy the boundary conditions B 



( 1 > 
±— 8 
l 2 ) 



= B . (b) ForS<<A. L , 



cosh — = 1 + — 
X 2 



( _ A 



, 8 . 1 
cosh — = 1 + — 
IX 2 



vVy 

r 8 ^ 



+ . 



\2Xj 



+ .. 



therefore B(x) = B a -B a (l/8Y)(8 2 -4x 2 ). 



2a. From (4), dF s = -MdB atT = 0. From Problem lb. 



M(x) = — — rB -(S 2 -4x 2 ), 

V ’ 4tt 8X 2 a 1 ’ 



whence 



F s (x,B a )-F s (0) 



1 

64t^ 2 



(S 2 -4x 2 )B ;i 2 . 



b. The average involves 



whence 



c. Let us set 



U 1, 3 ’ 



AF) = — B a 2 

96n 



f 8^ 2 

\Xj 



, for 8 « X. 



1 

96k 



B: 






— b; 

Sn 
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where B a f is the critical field for the film and B ac is the bulk critical field. Then 



B f = Jl2~B 

a f § 



3a. (CGS) curl H = ^ + - — 
c c <9t 



4n 

c 



a 0 E 



c 3 1 8E 

-A H . 

4nX~ ) c <9t 



curl curl H = -V 2 H = 



4n 

c 



a 0 curl E 

V 



c 3 1 dcurl E 

-B + 

4nh ) c 5t 



Now in CGS in nonmagnetic material B and H are identical. We use this and we use the 
Maxwell equation 



curlE = 



18 B 

c at 



to obtain 



V 2 B 



_Lr— — 

x 2 c 2 at 2 



4rar 0 5B 
c 2 at" 



If B ~ e 1 (k • r - cot) , then 



2 1 co" 47ua n co 



“ k "~FV + - 



= 0 . 



Q.E.D. 



b. 



1 

I 2 



2 2 , . 

CO CO 2 1 4rar 0 co 

— 1 — »—r\ also, co“ « 47ra 0 co and — » 

c" c” X c~ 



Thus the normal electrons play no role in the dispersion relation in the low frequency 
range. 

4. The magnetic influence of the core may be described by adding the two-dimensional 
delta function ® 0 8(p), where cj)o is the flux quantum. If the magnetic field is parallel to 
the z axis and div B = 0, then 



B -X 2 V 2 B = ® 0 8(p), 



or 
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B = -® 0 8(p). 



2 ( d 2 B | 1 gB" 

Up 2 + p ap y 

This equation has the solution B(p) = (® 0 /27rA, 2 )K 0 (p/A,), where K 0 is a hyperbolic 
Bessel function infinite at the origin and zero at infinity: 

(p«A,)B(p) = (® 0 /27iA, 2 )^n(A,/p); 

(p » A,)B(p)- (® 0 /27iA, 2 )(7iA./2p) 1/- exp(-p/^). 

The total flux is the flux quantum: 

oo oo 

2ttJ dp p B (p) = ® 0 Jdx x K 0 (x) = ® 0 . 

0 0 

5. It is a standard result of mechanics that E = -grad(p-c~'<9 A/<3t. If grad (p = 0, when 
we differentiate the London equation we obtain 0j/5t = (c 2 / 4nX L ~ j E. Now j = nqv and 

dj/dt = nq<3v/<9t = (nq 2 /m)E. Compare the two equations for cj/St to find 
c 2 / Aiik^ = nq 2 /m. 

Handbook of mathematical functions , U.S. National Bureau of Standards AMS 55, sec. 
9.6. 

6. Let x be the coordinate in the plane of the junction and normal to B, with 
-w/2 < x < w/2. The flux through a rectangle of width 2x and thickness T is 2xTB = 
(j) (x). The current through two elements at x and -x, each of width dx is 

dJ = (J 0 /w) cos[e® (x)/lic]dx = (J 0 /w)cos(2xTeB/Mc)dx, 
and the total current is 

/ W f 2 , sin(wTBe/Mc) 

J = (j„/w) | cos(xTeB/Hc)dx = J 0 . 

7a. For a sphere H (inside) = B ;i -4 tt:M/ 3; for the Meissner effect H (inside) = -4 tiM, 
whence B a =-8 ttM/3. 

b. The external field due to the sphere is that of a dipole of moment p = MV, when V is 
the volume. In the equatorial plane at the surface of the sphere the field of the sphere is 
-p/ a 3 = -4 tiM/ 3 = B a /2. The total field in this position is 3B a /2. 
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Chapter 11 



1. From Eq. (10), 



X = 




< r" 



> . 



i 00 _ / 

2 1 f 2 _2r / a o 2 

Here < r > = ^ • 4 ti r dr • e = 3a 0 . 

na 0 J 0 

23 —1 

The numerical result follows on using N = 6.02 x 10 mol . 

2a. Eu ++ has a half-filled f shell. Thus S = 7 x 1/2 = 7/2. The orbitals niL = 3, 2, 1, 0, -1, 
-2, -3 have one spin orientation filled, so that L = Zm L = 0. Also J = L + S = 7/2. Hence 
the ground state is 8 S 7 / 2 . 

b. Yb +++ has 13 electrons in the f shell, leaving one hole in the otherwise filled shell. 
Thus L = 3, S = 1/2, J = 7/2 — we add S to L if the shell is more than half-filled. The 
ground state symbol is 2 F 7 / 2 . 

c. Tb +++ has 8 f electrons, or one more than Eu ++ . Thus L = 3; S = 7/2 - 1/2 = 3; and J = 
6. The ground state is 7 F6. 

3a. The relative occupancy probabilities are 

e fA+ ‘ lB, / kT (Here A stands for k B A) 

e~ A/kT 

e -(A-nB)/kT 



1 

The average magnetic moment is 



< p > = p 



e -(A-MB)/kT _ e -(A+nB)/kT 



where Z = i + e ^ B )/ kT + e ^ kT +e - (A ^ B) / kT . 



b. At high temperatures e vkT — > 1 and 
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< p > — > p • 



pB + > 






kT -J 




l kT J 



4 



jli 2 B Np 2 

— • y — > — - — 

2kT ’ 2kT 



c. The energy levels as a function of field are: 




If the field is applied to take the system from a to b we 
increase the entropy of the spin system from « 0 to « 
N log 2. If the magnetization is carried out constant 
total entropy, it is necessary that the lattice entropy be 
reduced, which means the temperature -l . 



4a. Z = l + e~ A/T ; 

k.Ae"^ k„A 
l + e-vr e vr + i 



c = 



f dE^ 



v5T A 



= k B A 



4e A/T 

T 2 




b. For A/T « 1, e A/T -> 1 and C -> -k. 



i -T 
4 



5a. If the concentration in the spin-up band is N + = 1/2 N (1 + Q, the kinetic energy of all 
the electrons in that band is 



| N *£(3, t 3N *f = E 0 ( i +i ;f, 

and the magnetic energy is -N + p B = - 1/2 N(1 + Q p B. 
b. Now E tM = E 0 |(l + ^) 5 ' 3 +(1-^) 5 ,? }-N^iB; 
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^- = |e 0 {(i+<;) !,! -(i-(;) 2,3 }-%b = o 



■-C 

3 






9NpB _ 3pB 
20E 0 2s P 



M = NpC = ^*-B. Q.E.D. 

2s F 

6a. The number of pairs of electrons with parallel spin up is 

I(N*) J =i N ’( l + <;)\ 

so that the exchange energy among the up spins is 

4 VNl ( 1+ tf ; 

and among the down spins the exchange energy is 



■rvN 2 (i-Q 2 . 



b. Using these results and those from Prob. 5 we have E tot = E 0 |(l + ^) 5/3 + (l-<^) 5/J 
- - VN 2 (l + C 2 ) 2 - NCpB. Thus (for £ « 1 ) 



E 0 C--VN 2 C-NpB = 0 ; 
5C 9 0 2 



C = 



NpB 



— E 0 - — VN 2 
9 0 2 

NpB 

2Ns P _1 VN 2 



and 



m = Np^; = ■ 



3Np 2 



2s f --VN 
F 2 



-B 
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c. For B = 0 and C, = 0. 



d 2 E tot 20 E 

dCr 9 1 



- — VN 2 < 0 if V>— ^ 
2 9 N 2 



4 e F 
3~N 



7a. The Boltzmann factor gives directly, with x = keT 

e A A _ e - A A 

U = -A— 7 ji~ = -A tanh A/x; 

e A/T +e‘ A/T 7 

C = k B dU/dx = k B (A/k B T)“sech 2 (A/k B T), 



'y 

because d tanh x/dx = seclrx. 



b. The probability P(A) dA that the upper energy level lies between A and A + dA, 
referred to the midpoint as the zero of energy, is P(A) dA = (dA) / Aq. Thus, from (a), 



< U >= — J dA (A/A 0 ) tanh A/x, 

o 

Ao 

< C > = k B J dA (a 2 /a o x 2 j sech 2 (A/t) 

o 



x o 

= (k B x/A 0 )J dx x 2 sech 2 x , 

o 



where x = A/x. The integrand is dominated by contributions from 0 < A < x, because sech 
x decreases exponentially for large values of x. Thus 

oo 

< C > - (k B x/A 0 ) J dxx 2 seclrx. 

o 



_ <p> e IlB/T -e !lB/T 2 sinh x 

o. — — 

p 1 + 2 cosh x 1 + 2 cosh x 
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Chapter 12 



1. We have S e+§ Thus 



dS e X _ f 2JS 
dt i h 



6-T^\S y 



( 2JS^ 



v M j 



dS„ 



^ 2JS^ 



dt 



\ n j 



6-2(cos k x a + cos k y a + cos k z a) S p y ; 
6-2(cos k x a + cos k y a + cos k z a) 



S p • 



These equations have a solution with time-dependence ~ exp(-icot) if 
co = (2JS/M)(6-2 cos k x a-2 cos k y a-2 cos k z a). 

2. U = Z n k Mco k =hjdcoD(co)co<n(co)>. If co = Ak 2 , then dco/dk = 2Ak = 2 >/a co 1,/2 , 



and 



D(®) = 



4n co 1 



1 co 1/2 



87c 3 A 2VAco 1/2 47t 2 A 



3/2 



Then 



U= M 3/ J dco co 3/2 
4ttA 3/2 J 



1 



e pKco -r 



At low temps, 



r 1 r , x 3/2 _ 1 

- 1 : e“-l I hi;)” i“" 



(W 





r 


[Ad 


Uj 


s 

zeta 


U’ J 



gamma 
function v ' function 



[See Dwight 860.39] 

U = 0.45 (k B T) 5/2 / n 2 A 3/2 M 3/2 
C = dU/dT = 0. 1 13 k B (k B T/MA) 3/2 . 
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3 . 

M a T = C(B-pM B -sM a ) (B = applied field) 
M B T = C(B-sM B -pM A ) 



Non-trivial solution for B = 0 if 



T + sC pC 

pC T + sC 



0; T c = C(p-s) 



Now find % = (M a + M b )/B at T>T C : 



MT = 2CH-CM(s + p); % = ^ 

V x + c ^ + £ ) 

••• 0/T c = (p + e)/(p-e) • 

4. The terms in U e< + U c + U K which involve e xx are 



\_ 

2 



C u e xx - 




+ e„ 



) + Bia, 2 e xx . 



Take d/de xx : 



C n e xx + C 12 (eyy +e zz ) + B,a | 2 =0, for minimum. 



Further: 

C„e yy +C l2 (e xx +e zz ) + B,a 2 2 =0 . 
C 1 ie zz +C 12 (e xx +e yy ) + B 1 a 3 =0. 



Solve this set of equations for e xx : 

e =B C 12 -« Z : (C„ +2C 12 ) 
” 1 (C,, -C 12 )(C„ + 2C, 2 ) 

Similarly for e yy , e zz , and by identical method for e xy , etc. 
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5a. 

U (0) = K sin 2 0 - B a M s cos 0 

- Kcp 2 - B a M s (p 2 , for 0 = n + cp 

and expanding about small (p . 

For minimum near (p = 0 we need K>^-B a M s . Thus at B a = 2K/M S the magnetization 
reverses direction (we assume the magnetization reverses uniformly!). 

b. If we neglect the magnetic energy of the bidomain particle, the energies of the single 
and bidomain particles will be roughly equal when 

M s 2 d 3 ~a w d 2 ; or d c «a w /M s 2 . 

For Co the wall energy will be higher than for iron roughly in the ratio of the (anisotropy 
constant Kj) 1/2 , or VlO. Thus a w « 3 ergs/cm" . For Co, M s = 1400 (at room 

temperature), so M s 2 «2xl0 6 erg/cm 3 . We have dc«3/2xl0 6 «1.10~ 6 cm, or -100 A, 

as the critical size. The estimate is very rough (the wall thickness is d c ; the mag. en. is 
handled crudely). 

6 . Use the units of Eq. (9), and expand 

tanh — = — + [Dwight 657.3] 
t t 3 t 3 



_ _ , m m 

Then (9) becomes m- -+ ••• ; 

t 3t 3 



3(t 3 -t 2 ): 



2 2 ^.2 

m ; m 



3t 2 (l-t) 



but 1 - 1 is proportional to T c - T, so that moo^T, -T for T just below T c 



7. The maximum demagnetization field in a Neel wall is -4 7 tM s , and the maximum self- 
energy density is — ( 47 tM s )M s . In a wall of thickness Na, where a is the lattice constant, 



the demagnetization contribution to the surface energy is a dcrnag « 27iM s 2 Na . The total 
wall energy, exchange + demag, is a w « [ 71 2 JS 2 /Na 3 ) + ( 27iM s 2 Naj, by use of (56). The 
minimum is at 
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and is given by 



5g w /5N = 0 = — tc^ JS 2 /l\fa 2 + 27tM s 2 a, or 

( i , y / 2 

N = -jiJS 2 /M s V , 

V 2 J 



g w ~7iM s S(27rJ/a) 1/2 * (I0)(l0 3 )(l0 4 /l0 8 )' /2 « 10 erg/cm 2 , 

which is larger than (59) for iron. (According to Table 8.1 of the book by R. M. White 
and T. H. Geballe, the Bloch wall thickness in Permalloy is 16 times that in iron; this 
large value of 8 favors the changeover to Neel walls in thin films.) 

8. (a) Consider the resistance of the up and down spins separately. 

Magnetizations parallel: 

R n (up) = a p (LI A) + a; 1 (LI A) = 2a]' (LI A) 

(down) = cr fl 1 (L/ A) + a]' (LI A) = 2a] 1 (LI A) 

These resistances add in parallel: 

R n = R^ (down) R n (up) I[R n (down) + R n (up)] = 2 (L/A)/(a a +cr p ) 
Magnetizations antiparallel: 

R n (up) = a p 1 (L / A) + cr; 1 (L / A) 

R^^(clown) = a a '(LI A) + a p '(LI A) = R n (up) 

These (equal) resistances add in parallel : 

&U = R n («/>)■ 1 2 = (Li / A)(a; 1 + a p l ) / 2 
The GMRR is then: 

GMRR = R n IR n - 1 = (a; 1 + a p )(a u +a p )l 4-1 
= (a a /a p +a p /a a -2)/4 

(b) For the 'l -l magnetization configuration, an electron of a given spin direction must 
always go through a region where it is antiparallel to the magnetization. If a a — > 0, then 
the conductance is blocked and the resistance is infinite. 
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Chapter 13 



1. Consider a coil which when empty has resistance Ro and inductance Lo. The 
impedance is Zo = Ro - icoLo. When the coil is filled with material of permeability 
p = 1 + 47 1 % the impedance is Z = R 0 -icoL 0 (l + 47tx) = R 0 -icoL 0 (l + 47tx , + 47U %") , or 

Z = R 0 + 47t cox ff L 0 ~ icoL 0 (l + 4 tcx') • 

R 



dF dF „ dx 

2a. — = — -x + F y — + • 
dt dt ~ dt 



dF^ 




dx „ dy dz 


— — 


+ 




dt ) 


R 


y dt y dt 2 dt J 



Now 



dx , M dy , v dz ( M 
dt V ~ dt dt V ~ 



F — 4 = Q x F . 

dt ~ " 



, dM -- ~ 

b. — - = yMxB 
dt ~ ~ 



C dM ' 



V dt y R 



+ Q x M = yM x B 



^ dM^ 



V dt y R 



= yMx 



f 

Y J 



c. With Q = -yB 0 z we have 



dM 

dt 



/ R 



yM x B t x , 



so that M precesses about x with a frequency co = yBi. The time ti/ 2 to give ti /200 = 7T is 
1 1/2 = 7t/yB 1 . 

d. The field B, rotates in the xy plane with frequency Q = yB„. 
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3a. < B ; > = 



/ \ 2 

' a ' 

vNy 



ZZ< Ij Z I k z > , where for I = we have< Ij Z I k z > = ^-S jk . Thus 



b. 



< B; > = 



^ a V 1 



V Ny 



— Z5 jk = — — 
4 jk Jk 4N 



< B. > = 



f \ 4 
' a ' 

vNy 



z <i;i k %v>- 

jk&n J 



Now 



<I i 'Ik I I,‘I.‘>=^[S J1 5„8,„ + 8 Jl 8,„ 

to 

+8 j< 8 km +8 jm 8 M ], and 



< B t 4 > = 


fa} 


1 2 

— [N + 3N 2 ]^ 






In J 


16 


In, 



3N 2 

leT 



4. For small 9, we have U K - K9 2 . Now the magnetic energy density U M = -BM cos 9 - 

1 2 

-BM + — BM9“ , so that with proper choice of the zero of energy the anisotropy energy is 
equivalent to a field 

B a = 2K/M 

along the z axis. This is valid for 9 « 1 . For a sphere the demagnetizing field is parallel 
to M and exerts no torque on the spin system. Thus Bo + Ba is the effective field. 

5. We may rewrite (48) with appropriate changes in M, and with B an i sotrop y = 9. Thus 

-icoM/ =-iy A (M A + X|M B | + M B + X|M A |) ; 

-icoM B + =iy B (M £ 

The secular equation is 

Y A Z|M B | -CO 
-Yb^| m b| 

or 



^|m a | + m a + x|m b |) 



Ya^|m a | 

Yb^|M a |-co 



= 9 
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CO 2 - co(y A ^ |M B I - Yb^ | m a I) = 0 . 
One root is coo = 0; this is the uniform mode. The other root is 

®o = Mya I m b|-Yb| m a|)=°; 

this is the exchange mode. 
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Chapter 14 



1. E x0 = — — = kA sin kx e kz , and at the boundary this is equal to E xi . The normal 

dx 

component of D at the boundary, but outside the medium, is s(co)kA cos kx, where for a 
plasma s(co) = 1 - co p 2 /co 2 . The boundary condition is -kA cos kx = s(co)kA cos kx, or 
s(co) = -1, or co p 2 = 2co 2 . This frequency co = co p /V2 is that of a surface plasmon. 



2. A solution below the interface is of the form cp(-) = A cos kx e kz , and above the 
interface cp(+) = Acos kx c kz , just as for Prob. (1). The condition that the normal 
component of D be continuous across the interface reduces to Si(co) = -82(03), or 



1 — 



03 



03 



= -l + - 



03 



P 2 



03" 



, so that 03 2 = ~(®pi 



+ ®p2) 



3. (a) The equation of motion of the electrons is 

-G3 2 x e = -(e/m e )E x +icoo3 e y e ; -03 2 y e = -(e/m e )E y -io303 e x e . For the holes, 
-03 2 x h =(e/m h )E x +io303 h y h ; -03 2 y h = (e/m h )E y -icoo3 h x h . 



The result follows on forming £, e = x e + iy e and £ h = Xh + iyh- (b) Expand 

(o3 e + 03) 1 — co e 1 (1 — 03 / G3 e ) and (o3 h —03) — 03 h 1 (1 + 03 / co h ) . In this approximation 

(S h - ) / E + - (c/B)(co h 1 + co e 1 ) = (c 2 / eB 2 ) (m h + m e ) . 



4. From the solution to Problem 3 we have P + = pe E + / m h co h 03 , where we have dropped 
a term in co in comparison with cohco. The dielectric constant 

s(co) = 1 + 47 tP + / E + -47rpe 2 / m h co h co , and the dispersion relation s(co)co 2 = c 2 k 2 becomes 
47rpe 2 co/(eB/c) = c 2 k 2 . Numerically, co * [(10 3 )(3xl0 10 )/(10)(3xl0 22 )(5xl0 10 )] * 0.2 s 1 . 

It is true that cor will be «1 for any reasonable relaxation time, but co c x > 1 can be 
shown to be the applicable criterion for helicon resonance. 

5. md 2 r/dt 2 = -mco 2 r = -eE = 4a:eP/3 = -47ine 2 r/3 . Thus co 2 = 47me 2 /3m. 

6. md 2 r / dt 2 = -mco 2 r = -(e/c) (v x Bz) - mco o 2 r, where co 0 2 = 47ine 2 /3m, from the 
solution to A. Thus, with co c = eB/mc, 

-co 2 x = icoco c y - co o 2 x ; 

-co 2 y = -icoco c x - co 0 2 y . 
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2 2 ? 2 

Form £, = x + iy ; then -go £, - coco c £, + co 0 “£, = 0, or ok + coco c - co Q = 0, a quadratic 
equation for co. 

7. Eq. (53) becomes c 2 K 2 E = co 2 [s(oo)E + 4 tiPJ, where P is the ionic contribution to the 
polarization. Then (55) becomes 

co 2 s(oo)-c 2 K 2 4tt:co 2 
Nq 2 /M co 2 -co T 2 

or 

co 4 s 2 (oo)-co 2 c 2 K 2 +s(co)co T 2 +47iNq 2 /Mj + c 2 K 2 co T 2 =0. 

One root at K = 0 is co 2 = co T 2 + 47tNq 2 /s(oo)M. For the root at low co and K we neglect 
terms in co and in oTK". Then 

co 2 =c 2 K 2 co t 2 /[s(oo)co t 2 +47rNq 2 /M] 

= c 2 K 2 /[s(oo) +47rNq 2 /Mco T 2 ] = c 2 K 2 /s(0) , 

where s (0) is given by (58) with co = 0. 

8(a). a = ne 2 i/m = (co p 2 /4tt)t = 0.73 x 10 15 s' 1 = 800 (Q cm)- 1 

(b) co p 2 = 47tne 2 /m*; m* = 47ine 2 /co p “ = 4.2 xlO 27 g; m*/m = 4.7. 

9. The kinetic energy of a Fermi gas of N electrons in volume V is 

U = N(3/5)(M 2 /2m)(37r 2 N/V) 2/3 . Then dU/dV = - (2/3)U/V and d 2 U/dV 2 = (10/9)U/V 2 . 

The bulk modulus B = Vd 2 U/dV 2 = (10/9) U/V = (10/9) (3/5)n(mv F 2 /2) = nmv F 2 /3. 
The velocity of sound v = (B/p) 1,/2 , where the density p = n (m + M) - nM, whence 
v - (m/3M) 1/2 v F . 

10. The response is given, with p = 1/r, by 

m(d 2 x/dt 2 +pdx/dt + co p 2 xj = F(t) . 

The conductivity a does not enter this equation directly, although it may be written as a 
= co p 2 t/4ti. For order of magnitude, 
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a = (l/l0‘ 6 )(9xl0 n )-10 18 s 1 ; 
p = l/ x = v F /f = (l.6xl0 8 )/(4xl(T 6 )-0.4xl0 14 s " 1 ; 

co p = (47ine 2 /m) 1/2 - (lO x 10 23 x 23 x 10 20 /lO 27 )‘ /2 
- 1.5 x 10 16 s' . 



The homogeneous equation has a solution of the form x(t > 0) = Ae /J sin (cot + ((>), 



where co = 



i 1 / 2 



co„ 



+ (p/2) and X - p/2. To this we add the particular solution x 



-e/mco and find A and cj) to satisfy the initial conditions x(0) = 0 and x(0) = 0. 



11. The Laplacian V 2 cp = 0, whence 



d 2 f 

dz 2 



K 2 f = 0 . 



This has solutions 

f = Ae Kz for z < 0 
f = Ae~ K(z ~ d) for z > d 
f = B cosh K(z-d/2) for 0 < z < d . 



This solution assures that (p will be continuous across the boundaries if B = 
A/cosh(Kd/2). To arrange that the normal component of D is continuous, we need s(co) 
cxp/dz continuous, or s(co) = - tanh(Kd/2). 
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Chapter 15 



la. The displacement under this force is 

x ( t )=^£ a ( (0 ) e " cotdco - 

With co = cor + icoi, the integral is Ja(co)e K ' jR 'e fl>rt doj. This integral is zero for t < 0 

because we may then complete a contour with a semicircle in the upper half-plane, over 
which semicircle the integral vanishes. The integral over the entire contour is zero 
because a(co) is analytic in the upper half-plane. Therefore x(t) = 0 for t < 0. 

lb. We want 



x 



(•) 



1 f e~‘“ do 

O I 2 2 ’ 

27T^oCO 0 -co -Hop 



(A) 



which is called the retarded Green’s function of the problem. We can complete a contour 
integral by adding to x(t) the integral around an infinite semicircle in the upper half- 
plane. The complete contour integral vanishes because the integrand is analytic 
everywhere within the contour. But the integral over the infinite semicircle vanishes at t 
< 0, for then 



exp 



— i(co R + icOj ) ( — 1 1 ) = exp(-cOj |t|) exp(ico R t), 



which -A 0 as Icol — » <x>. Thus the integral in (A) must also vanish. For t > 0 we can 
evaluate x(t) by carrying out a Cauchy integral in the lower half-plane. The residues at 
the poles are 



±1 



(co 0 2 - jp 2 )' exp(-Tpt) exp +i(co 0 2 - jp 2 ) 



so that 



,(t) = (co 0 2 - IP 2 )' exp(-|pt)sin(co 0 2 - jp 2 )' 



2. In the limit co — » oo we have 



a 



'(«)-> -I f J/ 



co 



from (9), while from (11a) 
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J sa"(s)ds. 

0 

3. The reflected wave in vacuum may be written as 

-E y (reft) = B z (refl) = A' e~ i(kx+l0t) , 

where the sign of E y has been reversed relative to B z in order that the direction of energy 
flux (Poynting vector) be reversed in the reflected wave from that in the incident wave. 
For the transmitted wave in the dielectric medium we find 

E y (trans) = ck B z (transj/sco 
= s- 1/2 B z (trans) = A”e l(kx ’ tot) , 

2 2 2 

by use of the Maxwell equation c curl H = s5E/5t and the dispersion relation sco = c k 
for electromagnetic waves. 

The boundary conditions at the interface at x = 0 are that E y should be continuous: E y 
(inc) + E y (refl) = E y (trans), or A - A' = A". Also B z should be continuous, so that A + A' 
= s 1/2 A". We solve for the ratio A'/A to obtain s 1/2 (A - A') = A + A', whence 

A' _ 1 - s 1/2 
A “^7T 




and 



E(refl) _ A' _ s 1/2 -1 _ n + ik-1 
E(inc) A s 1/2 + 1 n + ik + 1 



The power reflectance is 



R(co) = r *r = 



f n-ik-1 2 ! 


( n + ik-l 2 ! 


v n-ik + l y 


v n + ik + l y 



(n-l) 2 +K 2 
(n + l) 2 + K 2 ’ 



4. (a) From (11) we have 



a 




ApfAhl 

r J 2 2 

n q s — co 



ds. 



In the limit co — » oo the denominator comes out of the integrand and we have 
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lim 

a 

co — » oo 



if 



(co) 



_ 2 _ 

7CC0 



J a'(s) ds. 

0 



(b) A superconductor has infinite conductivity at zero frequency and zero conductivity at 
frequencies up to co g at the energy gap. We can replace the lost portion of the integral 
(approximately a' n co g ) by a delta function a' n co g S(co) in a' s (co) at the origin. Then the KK 
relation above gives 



a " s ( co ) = ^: a 'n °Y 

7TC0 

(c) At very high frequencies the drift velocity of the conduction electrons satisfies the 
free electron equation of motion 

mdv/dt = -eE; -icomv = -eE, 



so that the current density is 

j = n (-e) v = -ine 2 E/mco 

2 

and coa" (co) = ne7m in this limit. Then use (a) to obtain the desired result. 



5. From (1 la) we have 



8 (co) — 1 



m J 0 s -co co g — co” 



22 1718 

6. n - K + 2inK = 1 + 4 tu Go/co. For normal metals at room temperature Go ~ 10 - 10 
sec -1 , so that in the infrared co<scg 0 . Thus n 2 -K 2 , so that R-l-2/n and 

n - ^(27tG 0 /co), whence R - 1 - ^(2co/7tg 0 ). (The units of g 0 are sec -1 in CGS.) 



7. The ground state of the line may be written v|/ g = A ] B 1 A 1 B 2 ...A N B N . Let the asterisk 
denote excited state; then if specific single atoms are excited the states are 
cpj =A,B 1 A 2 B 2 ...A j "B i ...A N B N ; 9j =A 1 B 1 A 2 B 2 ...A j B j \..A N B N . The hamiltonian acts 
thusly: 

7/cpj =8 A cp j +T l 0 j +T 2 0 H ; 

H 0j = s B 0 j + TjCpj + T 2 cpj + 1 . 
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An eigenstate for a single excitation will be of the form v|/ k = (acpj +(30j). We 

j 

form 



# V k = +aT,0 j +aT 2 e H 

j 

+ P £ b j + PTjtpj +PT 2 cpj +1]. 

= ^e« t -[(ae A +pT l +e-' l -pT,)< Pj 
j 

+ (aT,+p8 B +e ika aT 2 )0 j ] 

= E ¥k =Xe ijka [aE(p J +pE0 J ], 



This is satisfied if 

(s A -E)a + (T 1+ e- lka T 2 )p = 0; 
(T, + e ika T 2 )a + (s B -E)P = 0. 



The eigenvalues are the roots of 



s a -E T,+e- lka T 2 
T, + e ika T 2 s b -E 
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Chapter 16 



1. 



c X 3 , 3 3 

= eE; ex = r E = p; a = p/E = r = a H . 

r r 



4 n 4tc 

2. Ej = E 0 — — P = 0 inside a conducting sphere. Thus p = — a 3 P 

a = p/E 0 = a 3 . 



a 3 E 0 , and 



3. Because the normal component of D is continuous across a boundary, E air = sE die i, 
where E a ; r = 4 tiQ/A. with Q the charge on the boundary. The potential drop between the 

f O 



two plates is E air qd + E diel d = E air d 



q+- 

e; 



. For a plate of area A, the capacitance is 



C = 



A 



A ‘ 



47rd 



q+ 



e; 



It is useful to define an effective dielectric constant by 



1 1 

= — + q . 

£„ ff £ 



If s = oo, then s e ff = 1/q. We cannot have a higher effective dielectric constant than 1/q. 
For q = 1(T 3 , s e ff = 10 3 . 

4. The potential drop between the plates is Ei d + E 2 qd. The charge density 



Q _ D, _ sE, _ iq E 
A 4it 471 co 2 

by comparison of the way a and s enter the Maxwell equation for curl H. Thus 



(CGS) 



Ei + 



4ma 



SCO 



E, ; V = E,d 



4raa 
V sco 



+q 



Q = ^E 



CO 



2 ’ 



and thus C = — = 
V 



47rd 



_A 

1 icoq 
s 4rar _ 



and s 



eff 



=(i+q) T 



-(icosq/47ra) 
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4n 

F =F -P = F 

^int E 0 3 r E 0 



E. , = 



, 4n 

1h x 

3 



4n E 
— X int • 
3 4 



b-P = xE mt = 



, 4n 

1h x 

3 4 



6. E = 2Pi/a 3 . P 2 = aE = 2aPi/a 3 . This has solution pi = p 2 



0 if 2a = a 3 ; a = — a 3 . 



7 (a). One condition is, from (43), 

T(Tc-T ( ,)-|g 4 |P, 2 +g„P, 4 =0, 



The other condition is 



^y(T c -T„)P, 2 -t-|g 4 |P, 4 +i g# P, s = 0 



Thus 



-|g,|P, 2 +g 6 P. 4 =“|g4|P, 2 +jg 4 P, 4 ; 

— sP 2 =— k I ' p 2 =— 1®±1 

3 g 6 s 2 I§4|> s 4 gg • 



(b) From the first line of part (a), 



y(T c -T 0 ) 



3 |g 4 p 9 |g 4 

4 S 6 16 g 6 



T 2 
_3_g4_ 

16 g 6 



8. In an electric field the equilibrium condition becomes -E + y(T-T c )P + g 4 P 3 =0, 
where the term in gg is neglected for a second-order transition. Now let P = P s + AP . If 
we retain only linear terms in AP, then -E + y(T-T c )AP + g 4 3P s 2 AP = 0, with use of 
(40). Further, we can eliminate P s 2 because P s 2 = (y/g 4 )(T c -T) . Thus 
AP/E = l/2y (T c - T) . 
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9 a. |< — a — >| 
b. 1^ — 2a — ^1 




(na) 



Deforms to new stable structure of dimers, with lattice constant 2 x (former constant), 
c. 




10. The induced dipole moment on the atom at the origin is p = aE, where the electric 
field is that of all other dipoles: E = (2/a 3 )Zp n = (4p/a 3 )(Zn 3 );the sum is over 

positive integers. We assume all dipole moments equal to p. The self-consistency 
condition is that p = a(4p/a 3 ) (In 3 ), which has the solution p = 0 unless a > (a 3 /4) 

_o 

(1/In ). The value of the summation is 1.202; it is the zeta function q( 3). 
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Chapter 17 



1. (a) The interference condition for a linear lattice is a cos 9 = nA. The values of 9 that 
satisfy this condition each define a cone with axis parallel to the fiber axis and to the axis 
of the cylindrical film. Each cone intersects the film in a circle. When the film is flattened 
out, parallel lines result, (b) The intersection of a cone and a plane defines a conic 
section, here a hyperbola, (c) Let a, b be the primitive axes of a square lattice. The Laue 
equations (2.25) give a • Ak = 2 rap b • Ak = 2 7rr, where q, r are integers. Each equation 
defines a set of planes. The intersections of these planes gives a set of parallel lines, 
which play in diffraction from a two-dimensional structure the role played by reciprocal 
lattice points in diffraction from a three-dimensional structure. In the Ewald construction 
these lines intersect a sphere of radius k = 2 n/X in a set of points. In two dimensions any 
wavelength (below some maximum) will give points; in three dimensions only special 
values of X give points of intersection because one more Laue equation must be satisfied. 
The points correspond to the directions k' of the diffraction maxima. If the photographic 
plate is flat the diffraction pattern (2 dim.) will appear distorted. 



(ll)V 

i 

(To )'■ 



i-- -•r.-.d _ 

(ll) (0,1) (ID 




enter 

(beam direction) 



Points near the direction of the incident 
beam are shown. 



(d) The lattice of surface atoms in the 
(119) surface of an fee crystal is simple rectangular. The long side of the rectangle in 
crystal (real) space is a short side in the reciprocal lattice. This explains the 99° rotation 
between (21a) and (21b). 



2. With the trial function x exp (-ax), the normalization integral is 
| dx x 2 exp(-2ax) = l/4a 3 . The kinetic energy operator applied to the trial function 
gives 



-(li 2 /2m)d 2 u/dx 2 = -(h 2 /2m)(a 2 x - 2a)exp(-ax) 

while Vu = cEx“ exp (-ax). The definite integrals that are needed have the form 
| dx x 11 exp(-ax) = n!/a n+1 . The expectation value of the energy is < s >= (h 2 /2m)a 2 
+ (3eE/2a), which has an extremum with respect to the range parameter a when 
d < s >/da = (h 2 /2m)2a -3eE/2a 2 = 9, or a 3 =3eEm/2M 2 . The value of <s> is a 
minimum at this value of a, so that 

< £ > mm = (H 2 /2m)(3eEm/h 2 ) 2/3 +(3eE/2)(2h 2 /3eEm)' /3 
= (h 2 /2m )' /3 (3eE/2) 2/3 (2~ 2/3 + 2 1/3 ), 
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where the last factor has the value 1.89 .... The Airy function is treated in Sec. 10.4 of 
the NBS Handbook of mathematical functions. 



3. (a) D(s) = 



clN dk 
dk de 



2 d(7ik 2 ) m 
(2n I L) 2 dk Wk 



m 

~7dC- 



A 



2 

where A = L". 



Note: There are two flaws in the answer quoted in the text. First, the area A is 

missing, meaning the quoted answer is a density per unit area. This should not be a major 
issue. Second, the h should be replaced by fi . 



(b) N = 



2 

(2n/L) 2 



■nkj 



=> n s = N / A = kp !2n 



(c) R s = 



L m 
W n s e 2 r 



where n s is the 2D sheet density. For a square sample, W=L, so: 



R = 



2 n m 

Ti T~ 
k F e t 



and using hk F / m = v F : 

_ 2 nh _ h 1 

k F v F e 2 x e 2 k F t 
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CHAPTER 18 



1. Carbon nanotube band structure. 

(a) b ; - a ; - => b,=(-^^), b 2 =&,%). 

(b) The angle between K and biis 30° ; A right triangle is formed in the first BZ with two 

sides of length K and bill. Now bj = , so: 

K = (b,l 2)/cos(30°)= 47r/3a . 

(c) Quantization of k along x: k x (na)=2nj= k x =l7ij/na. 

Assume n = 3 i, where i is an integer. Then: k x = K(j/2i). For j = 2i, k x =K. Then 

AK = Lj and there is a massless subband. 
y J 

(d) For n = 10, k x =2n j/lOa =K(3j/20). The closest k comes to K is for j = 7, where 
Ak x = KZ20. Then: 

s n = 2hv F (4nl3a)/10 = 1.8 eV. 

The next closest is for j = 6, where Ak x = Kl 10, twice the previous one. Therefore: 

S '. 22 = 2S\\ . 

(e) For the lowest subband: |Ak|‘ = (K / 20) 2 + k 2 y , so: 

= [(HK / 20v F )v 2 F ] 2 +(hk y v F ) 2 



This is of the desired form, with m * = fiK / 20v f . 
m * / m = tiK / 20 mv F = 0.12. 

2. Filling subbands 

*2 2 *2 2 

s(n ,n ) = ~(n 2 + n 2 ) => States are filled up to £-(2,2)= — (8) 

1 y 2mW 2 x y 2 mW 2 



(1,1) subband: 



h 2 k 2 j 



2m 2mW 



h 2 7T 2 n , ^6 n 

r (8-2) => * u = — 



W 



2 , 2V6 

=> n \,\ ~ * 1,1 — u/ 
n Vv 



uu ^ t2k l\ h 2 ft 2 to . , 43n 2, 2V3 

(2,1) subband: — = - — — r (8-5) => k 2l = => n 2 l = — k 21 = - 



2m 2mW~ 



W 



n 



W 
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(2,1) subband: same. 



2V6 4V3 

n = 1 

W W 



5.9 x 10 8 /m. 



3. Breit-Wigner form of a transmission resonance 

(a) cos(8(p) = 1 - 8(p 2 / 2 ; I r. t 1=^/1- It, l 2 sl -{It, I 2 -|l( I 4 

The denominator of (29) is then: 

l + (l-lt, l 2 )(l- 1 r 2 I 2 ) -2(1 -{It, I 2 -jl t, l 4 )(l-{lt 2 1 2 -|lt 2 l 4 )(l-|^ 2 ) 
= {(It, I 4 + 1 1 2 l 4 ) + {lt, l 2 lt 2 1 2 +^ 2 = {(lt, I 2 + lt 2 1 2 ) 2 +8<p 2 
7 „. 4 1 1, I 2 1 1 2 I 2 

(It, I 2 +lt 2 1 2 ) 2 + 48(p 2 ' 



(b) 8cp = 2LSk and 8k / Ss = Ak / As = (a / L) / As . Combining: 
8(p = (2L)(ir / L)8s / => 8(p!2n = 8s / 



(c) Combining: 



4 It, l 2 l t 2 I 2 {As! In) 1 
(A£-/2tt) 2 (I t, I 2 + 1 1 2 1 2 ) 2 + 4c& 2 



which is (33). 



4. Barriers in series and Ohm’s law 

(a) 

1 _ 1- 1 r\ l 2 l r 2 1 2 _ t | 1- 1 r, 1 2 1 r 2 1 2 — 1 1, l 2 lt 2 l 2 
3 ~ It, l 2 lt 2 l 2 ~ + It, l 2 lt 2 l 2 

_ l | 1 — (I r 2 1 2 + 1 1 2 1 2 )+ 1 1, l 2 l r 2 I 2 +1 r, l 2 l t 2 I 2 
1 1, l 2 l t 2 1 2 



n 1 - (1- 1 1, I 2 ) I r 2 I 2 -(1- I r, I 2 ) 1 1 2 I 2 
It, l 2 lt 2 l 2 



= 1 + 



I A I 2 , I t\ I 



lt 2 I- 



+ - 



It, I 



which gives (36). 



(b) 



n, D e 2 r 
m 



2k F e 2 r 

mn 



and 




m 



=> 



2v F e 2 r _ 2e 2 (2v F r) 

fl7T h 



But: £ b =v f t b =2v f t 



=> 



<J, 



ID 



2e% 

h 



5. Energies of a spherical quantum dot 



(a) J E • da = Q enc , / ss o => E = q / 4 nss o r 2 



R+d 



v=l 



qdr 



9 



Anssr~ 4 nss 



1 



1 



0 v 



R R + d 



Integrating from inner to outer shell: 

q d 
4 nss o R( R + d) 
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_ q . R(R + d) tt e ~ e ~ 

C = — = 4 7Z££„ and therefore U = — = 



C 4 kss 0 R(R + d) 



(b) For d « R , 



„ . R 2 A 

C = 4 nee n — = ss — . 
° d ° d 



(c) For d» R , U = 



Akss.R 



. Also £,, n = ■ 



2m R 2 



U e 1 2m R 1 
£ 00 4-7T££ 0 R H 2 7T 2 



U _ e 2m R _ e 2m R 2 R 
£ 00 An££ 0 R fl 2 7V 2 4 n££ a ffTT 2 K 2 a B 



6 . Thermal properties in ID 

. . _ 2K 1 L 



(a) D(a>) = 



2 n!Lv m 



w r d(oD{o))fi<x> ~ fiL 


i oxl co tiL 


(k,T ' | 


2a 

1 1 


j- xdx 


' exp(/z<z>/ k B T)-\ 7W * 


\QXp(ha)l k B T)-\ TJX 


l h J 


J 

( 


' exp(.v) - 1 



Obtaining value from table of integrals: 

_ Lk\T 2 n 2 = n 2 Lk 2 B T 2 

fuw 6 3 hv 

^ xjj ,-rr I 2 7l 2 Lk 2 B T 

C v =dU tol /dT\ = 

3 hv 

(b) The heat flow to the right out of reservoir 1 is given by: 



J R 



D r (co) 



dcofico 



■1 2 ni 



H~s f k B T x V n 2 _ n 2 k 2 B T 2 



L exp(fta>/ k B T x )-l 2n ^ < 
and similarly for J/, The difference is: 



txpihcol k B T x ) -1 2n\ h ) 6 6 h 



J j ^i3( r 2 _ r 2) 
R L 6 h {I 2 ’ 



Let T x = T + AT , T 2 = T => (t 2 -T;)~ 2TAT for small AT. 



=> Jr-J l = 



71 ^ B ' S AT which gives (78). 
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Chapter 20 



1. U = nEi. 

The number of ways to pick n from N is N! / (N-n) !n! . The number of ways to put n 
into N' = n'! / (N' - n) !n! . 



N! N'! 

Entropy + j . 

l0g (N^ii)!n! ^ N1 ° gN ~( N ~ n ) l0g ( N ~ n )~ nl0gn 
l0g (N'-n)! n ! ~ N ' log N '-( N '- n ) lo g( N '- n )- nl °g n 

f dF ) au as .... . . 

— = T — = 0 in equilibrium; thus 

^an J T dn dn 

^ N-n , N'-n^ 

Ej = k B T log + log 

In n 



= k B T log 



(N-n)(N'-n) 



. For n « N, N' , 



Ej = k B T log (NN'/n 2 ) ; n = VNN 7 e _El 



-E,/2k B T 



2. From (2), n = Ne For sodium, N = 2.5 x 10 22 atoms/cm 3 . Thus 



n - ^2.5 x10 22 ^4x 10 18 ^«10 5 per cm 3 . 



IQ/: q 

3a. Mco = — ■ — eV. From Table 10.3 we have s (qo) = 2.25 for NaCl, whence 
s“ 4 

Mco = 2.0 eV. The observed value quoted in Table 18.2 is 2.7 eV. 



b. From Vol. I of Atomic energy levels . Circular 467 of the Nat. Bur. of Standards, 3p-3s 
= 16960 cm =2.1 eV. 
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Chapter 21 



1. In fee the nearest neighbors to a given atom lie in [110] and equivalent directions — 
thus these directions are lines of closest packing. In bcc the nearest neighbors lie in [111] 
and equivalent directions. There is a lively interest in the channeling of fast ions along 
open channel directions in crystals. 

2a. 



oooooooooo 
o ■■■ e— e— e~ o — o o o o o 

o o o o o ■ © - e~ e — o -e- < — remove these atoms 
oooooooooo 

b. 



remove — * 

these 

atoms 



ooooooooo 

disl. » o o o o o - o - o o - < — disl. 

ooooooooo 



3. When a dislocation moves a distance L through the crystal, the crystal advances by a 
distance b. The work aL“b may be expressed as 




FL • L 

force ^distance 
on moved 

disloc . 



Thus F = ob is the force/unit length on 
the dislocation. 
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Chapter 22 



1. Disordered state: S given by Eq. (2.48) with <f> as averaged over 3 Cu and 1 Au atom. 
Ordered state: 

S(hk(?) = f Au +f Cu (e~ III<k+,) +e~‘ (h+r) +e~ 1II<h+k *) Consider the following reflections: (x = 
refl. present) 



Indices 


Disordered 


Ordered 


100 


no 


X 


110 


no 


X 


111 


X 


X 


200 


X 


X 


210 


no 


X 


211 


no 


X 


220 


X 


X 


221 


no 


X 


222 


X 


X 



2. C = 5E/5T; now use E from Eq. (5): E = E 0 + 2NUP 2 . Thus C con f ig = 4NUP(T) OP/oT, 
and we recall that U is negative, as otherwise an ordered state does not occur. 

From Fig. 7b, 




in this region a treatment based on P alone gives 
Cconfig = 0. We should look at the short range order! 
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